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Abstract. In this article we study the pointwise decay properties of solu- 
tions to the wave equation on a class of nonstationary asymptotically flat 
backgrounds in three space dimensions. Under the assumption that uniform 
energy bounds and a weak form of local energy decay hold forward in time 
we establish a t~ 3 local uniform decay rate (Price's law [54]) for linear waves. 
As a corollary, we also prove Price's law for certain small perturbations of the 
Kerr metric. 

This result was previously established by the second author in [64 ) on sta- 
tionary backgrounds. The present work was motivated by the problem of 
nonlinear stability of the Kerr/Schwarzschild solutions for the vacuum Ein- 
stein equations, which seems to require a more robust approach to proving 
linear decay estimates. 



1. Introduction 

In this article we consider the question of pointwise decay for solutions to the 
wave equation on certain asymptotically flat backgrounds. Our interest in this prob- 
lem originates in general relativity, more precisely the wave equation on Schwarzs- 
child and Kerr backgrounds. There the expected local decay, heuristically derived 
by Price [54] in the Schwarzschild case, is i~ 3 . This conjecture was considered 
independently in two recent articles |26j and [64] . 

The work in [26] is devoted to the Schwarzschild space-time, where separation 
of variables can be used; in that context, very precise sharp local decay bounds arc 
established for each of the spherical modes. Precisely, it is shown that the A:-th 
mode leads to t~ 3 ~ 2k local decay. 

The work in |64j , on the other hand, applies to a large class of stationary asymp- 
totically flat space-times, and asserts that if local energy decay holds then Price's 
law holds. Further, sharp pointwise decay rates are established in the full forward 
light cone; these have the form t _1 (t — r)~ 2 . Local energy decay (described later 
in the paper) is known to hold in Schwarzschild and Kerr space-times. 

Both of the above results involve taking the Fourier transform in time and hence 
rely heavily on the stationarity assumption. The aim of this article is to prove 
the same result as in [64], namely that local energy decay implies Price's law, but 
without the stationarity assumption. The proof below is more robust than the one 
in [64] . and improves on the classical vector field method. As an application, in the 
last section of the paper we prove that local energy decay (and thus Price's law) 
holds for a class of nonstationary perturbations of the Kerr space-time. 
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Just as in [64] , this work is based on the idea that the local energy estimates con- 
tain all the important local information concerning the flow, and that only leaves 
the analysis near spatial infinity to be understood. In the context of asymptoti- 
cally flat metrics, this idea originates in earlier work [65| and |48| where it is proved 
that local energy decay implies Strichartz estimates in the asymptotically flat set- 
ting, first for the Schrodingcr equation and then for the wave equation. The same 
principle was exploited in }45) and |67) to prove Strichartz estimates for the wave 
equation on the Schwarzschild and then on the Kerr space-time. 

1.1. Notations and the regularity of the metric. We use (t = xq,x) for the 
coordinates in R 1+3 . We use Latin indices i, j = 1,2,3 for spatial summation and 
Greek indices a, j3 = 0, 1,2,3 for space-time summation. In R 3 we also use polar 
coordinates x = ru with u€S 2 . By (r) we denote a smooth radial function which 
agrees with r for large r and satisfies (r) > 2. We consider a partition of R 3 into 
the dyadic sets Ar = {(r) « R} for R > 1, with the obvious change for R = 1. 

To describe the regularity of the coefficients of the metric, we use the following 
sets of vector fields: 

t = di}, n = {x t d 3 - xjdi}, s = td t + xd m , 

namely the generators of translations, rotations and scaling. We set Z = {T, f2, S}. 
Then we define the classes S z (r k ) of functions in R + x R 3 by 

a g S z (r k ) \Z j a{t,x)\ < Cj (r) k , j > 0. 

By S z ad (r k ) we denote spherically symmetric functions in S z (r ). 

The estimates in this article apply to solutions for an inhomogeneous problem 
of the form 

(l.l) (n g + V)u = f, u(0) = u , am(o) = m 

associated to d'Alembertian D g corresponding to a Lorentzian metric g, a potential 
V, nonhomogeneous term / and compactly supported initial data uq, u\. For the 
metric g we consider two cases: 

Case A: g is a smooth Lorentzian metric in M + x R 3 , with the following prop- 
erties: 

(i) The level sets t = const are space-like. 

(ii) g is asymptotically flat in the following sense: 

g = m + g sr + gi r , 

where m stands for the Minkowski metric, gi r is a stationary long range spherically 
symmetric component, with ^^^(r -1 ) coefficients, of the form 

gw = girM{r)dt 2 + gi rttr {r)dtdr + gi r ,rr(r)dr 2 + gi rtUlul (r)r 2 duP 1 

and g sr is a short range component of the form 

g sr = g sr ,ttdt + 2g sr jidtdxi + gsr.ijdXidXj 

with S z (r~ 2 ) coefficients. 

We remark that these assumptions guarantee that dt is time-like near spatial 
infinity, but not necessarily in a compact set. This leads us to the second case we 
consider: 
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Case B: g is a smooth Lorentzian metric in an exterior domain KxR 3 \B(0,i?o) 
which satisfies (i),(h) above in its domain, and in addition 

(iii) the lateral boundary R x 8B(0,Rq) is outgoing space-like. 

This latter condition insures that the corresponding wave equation is well-posed 
forward in time. This assumption is satisfied in the case of the Schwarzschild and 
Kerr metrics (or small perturbations thereof) in suitable advanced time coordi- 
nates. There the parameter Ro is chosen so that < Ro < 2M in the case of the 
Schwarzschild metric, respectively r~ < Ro < r + in the case of Kerr (see [3U] for 
the definition of r^), so that the exterior of the Rq ball contains a neighbourhood 
of the event horizon. 

1.2. Normalized coordinates. Our decay results are expressed relative to the 
distance to the Minkowski null cone {t = \x\}. This can only be done provided 
that there is a null cone associated to the metric g which is within 0(1) of the 
Minkowski null cone. However, in general the long range component of the metric 
produces a logarithmic correction to the cone. This issue can be remedied via a 
change of coordinates which roughly corresponds to using Rcggc- Wheeler coordi- 
nates in Schwarzschild/Kerr near spatial infinity. See |64j . After a further conformal 
transformation^, see also |64j . the metric g is reduced to a normal form where 

gir = gu(r)r 2 du 2 , g^ £ Sf^r -1 ). 

In particular, we can replace O g by an operator of the form 

(1.2) P = U + Q 

where □ denotes the d'Alembertian in the Minkowski metric and the perturbation 
Q has the form 

(1.3) Q = sf^ + dag^dp + V, gZ P € S 2 ^ 2 ), gf e S^ ad {r^), V € S z {r^). 

We call these coordinates normal coordinates. All of the analysis in the paper is 
done in normal coordinates and with g in normal form. The full perturbation Q 
above has only short range effects. 

1.3. Uniform energy bounds. The Cauchy data at time t for the evolution (jTTTJ) 
is given by (u(i) , dtu(t)) . To measure it we use the Sobolev spaces H k , with the 
qualification that in Case A this means H k := H k (M. 3 ), while in Case B we use 
the obvious modification H k := H k (M. 3 \B(0, Rq)). Wc begin with the following 
definition: 

Definition 1.1. We say that the evolution (|1.1[) is forward bounded if the following 
estimates hold: 

(1.4) ||V«(ti)|| Hk <Cfc(||V«(to)IU* + ||/IUi([t 0l *i]iH*)), 0<t o <ti, k>0. 

It suffices to have this property when / = 0. Then the / term can be added in 
by the Duhamel formula. One case when the uniform forward bounds above are 
easy to establish is when d t is a Killing vector which is everywhere time-like and 
V is nonnegative and stationary. Otherwise, there is no general result, but various 
cases have been studied on a case by case basis. 

We remark that in the case of the Schwarzschild and Kerr space-times dt is not 
everywhere time-like, so the forward boundedness is not straightforward. However, 
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it is known to hold for Schwarzschild (see [24] and [45]) as well as for Kerr with 
small angular momentum (see [29], [23], [66]) and for a class of small stationary 
perturbations of Schwarzschild (see [23]). 

The forward boundedness is not explicitly used in what follows, but it is defined 
here since it is usually seen as a prerequisite for everything that follows. 

1.4. Local energy decay. A stronger property of the wave flow is local energy 
decay. We introduce the local energy norm LE 

\\u\\le = sup|j(r)~^w|| L 2 (RxAj?) 

(1-5) 

\W\\LE[t M] =sup||(r) 2 u\\ L 2 {[toM]xAR) , 

R 

its H 1 counterpart 

IMIlb 1 = l|Vu|| LB + || (r)- 1 ^!^ 

(1.6) 

IMUi5i[to,ti] = ll Vu llLB[t ,ti] + 11(0 U\\LE[toM> 

as well as the dual norm 

\\f\\LE> = Il( r >'/IU 2 (KXA R ) 

(1-7) * , 

ll/]|i£7-[t ,*i] = 2^ IK^) 2 /IU 2 ([to,ti]xA R )- 
R 

These definitions are specific to (1 + 3)-dimcnsional problems. Some appropriate 
modifications are needed in other dimensions, see for instance |48| . We also define 
similar norms for higher Sobolev regularity 

\\u\\ LEl , k = J2 \\d a u\\ LE i 

\a\<k 

WuWLE^toM = H^IU-EMto.ti] 
|a|<fc 

\\u\\LE0."[t , tl ] = H \\d a u\\ LE[toM] , 



\a\<k 



respectively 



\a\<k 

ll/l|i^*.*[to,*i] = X! ll 5 "/IUs«[t ,ti]- 
\a\<k 

In Case A above this leads to the following 

Definition 1.2. We say that the evolution (jl.l|) has the local energy decay property 
if the following estimate holds: 

(1.8) \\u\\ LE i, k[t0iOo) < c fc (||Vu(to)||jr* + ||/IUjs-.*[*o,«»))» *o>0, fc>0 
mix K 3 . 

The first local energy decay estimates for the wave equation were proved in the 
work of Morawetz [50] , [5T] , [52] ; estimates of this type are often called Morawetz 
estimates. There is by now an extensive literature devoted to this topic and its 
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applications; without being exhaustive we mention [59], [36], [56], [34], [35], [9], 
SI, S3, [H], [31]. 

The sharp form of the estimates as well as the notations above are from Mctcalfc- 
Tataru [48]; this paper also contains a proof of the local energy decay estimates 
for small (time dependent) long range perturbations of the Minkowski space-time 
and further references. See also [49], [I], [46] for time dependent perturbations, as 
well as, e.g., [8], [7], |57j for time indepedent, nontrapping perturbations. There is 
a related family of local energy estimates for the Schrodinger equation. See, e.g., 
the original works [T7], [55], [55] in this direction as well as [25J, [T5] in variable 
geometries. For notations and estimates most reminiscent to those used here, we 
refer the reader to [65] and [44] . 

In Case B an estimate such as ([1.8)1 cannot hold due to the existence of trapped 
rays, i.e. null geodesies confined to a compact spatial region. However a weaker form 
of the local energy decay may still hold if the trapped null geodesic are hyperbolic. 
This is the case for both the Schwarzschild metric and for the Kerr metric with 
angular momentum \a\ < M. To state such bounds we introduce a weaker version 
of the local energy decay norm 

\H\lei = ||(1 -x)V«||x,b + ||(r)- 1 «||iB 

NUfiilto.ti] = \\( l -X)^u\\ LElta.tr] + II (r) _1 «|| LE[t M] 

for some spatial cutoff function \ which is smooth and compactly supported. 
Heuristically, \ is chosen so that it equals 1 in a neighbourhood of the trapped 
set. We define as well a dual type norm 

ll/lks* = \\xf\\t»m + \\0—x)f\\LE* 

ll/IUE*t*o,ti] = llx/llL^to,*!]^ 1 + IK 1 - X)f\\hE*[t M}- 

As before we define the higher norms LE^; k respectively LE^ k . 

Definition 1.3. We say that the evolution (|1.1|) has the weak local energy decay 
property if the following estimate holds: 

(!-9) IMIl^W) ^ Ck(\\Vu(t )\\ H * + ||/IU«,* [t0iOo) ), k > 0, to > 
in either Ixl 3 or in the exterior domain case. 
Note that this implies in particular 

(1-10) IM| iS i,fc[i 0lOO ) < C k (\\S7u{t )\\ Hk +i + ||/||z,B.,fe+i[t ,oo), 

hence we can get rid of the loss at the trapped set by paying the price of one extra 
derivative. 

Two examples where weak local energy decay is known to hold are the Schwarz- 
schild space-time and the Kerr space-time with small angular momentum |a| <C M. 
In the Schwarzschild case the above form of the local energy bounds with k = 
was obtained in g5], following earlier results in SB, El, 0,0, ©, El' Thc 
number of derivatives lost in (|1.10D can be improved to any e > (see, for example, 
[3], S3), but that is not relevant for the problem at hand. 

In the case of the Kerr space-time with small angular momentum |a| <JC M the 
local energy estimates were first proved in Tataru-Tohancanu [66] , in a form which 
is compatible with Definition 11.31 Stronger bounds near the trapped set as well as 
Strichartz estimates are contained in the paper of Tohancanu [67] . For related work 
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we also refer the reader to [67], [20] and [2]. For large angular momentum \a\ < M 
a similar estimate was proved in |19| for axisymmetric solutions. 

The high frequency analysis of the dynamics near the hyperbolic trapped orbits 
is a very interesting related topic, but does not have much to do with the present 
article. For more information we refer the reader to [16], [14], [53], [70], [10] . 

One disadvantage of the bound (|1.10j) is that it is not very stable with respect to 
perturbations. To compensate for that, for the present result we need to introduce 
an additional local energy type bound: 

Definition 1.4. We say that the problem (|1.1|) satisfies stationary local energy 

decay bounds if on any time interval [to,ti] and k > we have 

(1.11) 

\\ u \\LE^"lt .t 1 ] \\Vu(t Q )\\ Hk + ||Vw(ti)|| H fc + ||/||L£*>*[t ,ti] + \\dtu\\ LEa ,k [toM] . 

Unlike the weak local energy decay, here there is no loss of derivatives. Instead, 
the price we pay is the local energy of d t u on the right. Heuristically, (|l.llj) can be 
viewed as a consequence of (|1.9|) whenever d t is timelike near the trapped set. In 
the stationary case, this can be thought of as a substitute of an elliptic estimate at 
zero frequency. 

While one can view the stationary local energy decay as a consequence of the 
local energy decay, it is in effect far more robust and easier to prove than the weak 
local energy decay provided that d t is timclikc near the trapped set. This allows one 
to completely sidestep all trapping related issues. This difference is quite apparent 
in the last section of the paper, where we separately establish both stationary local 
energy decay and weak local energy decay for small perturbations of Kerr. While 
the former requires merely smallness of the perturbation uniformly in time, the 
latter needs a much stronger t~ l ~ decay to Kerr. 

Because of the above considerations, for our first (and main) result in the theorem 
below we are only using as hypothesis the stationary local energy condition. 

1.5. The main result. Given a multiindex a we denote by u a = Z a u. By u< m 
we denote the collection of all u a with \a\ < m. We are now ready to state the 
main result of the paper: 

Theorem 1.5. Let g be a metric which satisfies the conditions (i), (ii) mix 
R 3 , or (i), (ii), (Hi) in M x M 3 \ B(0,R o ), and V belonging to S(r~ 3 ). Assume 
that the evolution (jl.ljl satisfies the stationary local energy bounds from Definition 
\l-4\ Suppose that in normalized coordinates the function u solves Ou = f and is 
supported in the forward cone C = {t > r — Ri} for some R\ > 0. Then the 
following estimate holds in normalized coordinates for large enough m: 

(1.12) \u(t,x)\< K {t){t 1 _ r)2 , \Vu(t,x)\< K ^—^ 

where 

K = \\u< m \\ LE i + \\t?f< m \\ LE ' + \\rt^Vf< m \\ LE '. 

If in addition the weak local energy bounds (|1.9[) hold then the same result is valid 
for all forward solutions to (jTTTJ) with data (uq,ui) and f supported inside the cone 
C and 

k = ||Vu(0)||ifm + f< m \\LE' + \\rt^f< m \\ LE «. 
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We remark that we actually prove a slightly stronger result, with k replaced by 
C4 in Lemma f3 . 2 H 

As an application of this result, in the last part of the paper we prove Price's 
Law for certain small, time-dependent perturbations of Kerr spacetimcs with small 
angular momentum (and V = 0). 

The problem of obtaining pointwise decay rates for linear and nonlinear wave 
equations has had a long history. Dispersive L 1 — > L°° estimates providing t" 1 
decay of 3 + 1 dimensional waves in the Minkowski setting have been known for a 
long time. 

The need for weighted decay inside the cone arose in John's proof [32] of the 
Strauss conjecture in 3 + 1 dimensions. Decay bounds for □ + V with V = 0(r~ 3 ), 
similar to those given by Price's heuristics, were obtained by Strauss-Tsutaya |60j 
and Szpak [62], [61]. See also Szpak-Bizoh-Chmaj-Rostworowski [63]. 

A more robust way of proving pointwise estimates via L 2 bounds and Sobolev 
inequalities was introduced in the work of Klainerman, who developed the so-called 
vector field method, see for instance |38j , |37j . This idea turned out to have a myriad 
of applications and played a key role in the Christodoulou-Klainerman [15| proof 
of the asymptotic stability of the Minkowski space time for the vacuum Einstein 
equations. 

In the context of the Schwarzschild space-time, Price was the first to heuristically 
compute the t~ 3 decay rate for linear waves. More precise heuristic computations 
were carried out later by Ching-Leung-Suen- Young [12], [13]. Following work of 
Wald [69], the first rigorous proof of the boundedness of the solutions to the wave 
equation was given in Kay- Wald [33] . 

Uniform pointwise i -1 decay estimates were obtained by Blue-Sterbenz [3] and 
also Dafermos-Rodnianski [24]; the bounds in the latter paper are stronger in that 
they extend uniformly up to the event horizon. A local t~? decay result was 
obtained by Luk [43]. These results are obtained using multiplier techniques, related 
to Klaincrman's vector field method; in particular the conformal multiplier plays a 
key role. 

Another venue which was explored was to use the spherical symmetry in order 
to produce an expansion into spherical modes and to study the corresponding ode. 
This was pursued by Kronthaler [40], [39], who in the latter article was able to 
establish the sharp Price's Law in the spherically symmetric case. A related analysis 
was carried out later by Donninger-Schlag-Soffer [27] for all the spherical modes; 
they establish a t~ 2 ~ 2k local decay for the fc-th spherical mode. Later the same 
authors obtain the sharp t~ 3 ~ 2k in [2"S] . 

Switching to Kerr, the first decay results there were obtained by Finster-Kamran- 
Smoller-Yau [28] . Later Dafermos-Rodnianski [20] and Anderson-Blue [2] were able 
to extend their Schwarzschild results to Kerr, obtaining almost aT 1 decay. This 
was improved to i -1 by Dafermos-Rodnianski [21] and to t~? by Luk |42j . 

Finally, the i~ 3 decay result (Price's Law) was proved by the second author 
in [M]; the result there applies for a large class of stationary, asymptotically flat 
space-times. In addition, in |64j the optimal decay is obtained in the full forward 
light cone. 
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2. Vector fields and local energy decay 

The primary goal of this section is to develop localized energy estimates when 
the vector fields Z are applied to the solution u. 

2.1. Vector fields: notations and definitions. For a triplet a = (i,j,k) of 
multi- indices i, j and nonnegative integer fc we denote |a| = |z| + 3|j| + 9k and 

u a = T i n ] S k u. 

On the family of such triplets we introduce the ordering induced by the ordering of 
integers. Namely, if u\ = (ii,ji,k\) and a 2 = (i2 3 J*2 j ^2) then we define 

oti < a 2 = \h\ < \ji\ < \h\, h < k 2 . 

We use < for the case when equality does not hold. For an integer m we denote 

0i\ < c*2 + m = 0i\ < a.2 + /3, |/3| < m. 

We also define 

^<m (^a)|a|<m 
U<f3 = {u a ) a <p. 

and the analogues for < instead of <. 

We now study the commutation properties of the vector fields with P. Denoting 
by Q sr the class of all operators of the form 

d a g af) d p + V, g^eS z (r- 2 ), V G S z (r~ 3 ), 

we see that Q defined in (|1.3[) consists of RVL 2 where R £ S z ad (r~ 3 ) plus an element 
of Q sr . 

We now record the commutators of P with vector fields. The commutator with 
T yields 

(2.13) [P,T] = Q, QeQ sr . 
The same applies with il, 

(2.14) [P,Q] = Q, QeQ sr . 

However, in the case of S we get an extra contribution arising from the long range 
part of P, 

(2.15) PS-{S + 2)P = Q + m 2 , QeQsr, i? G S z ad (r~ 3 ). 

Further commutations preserve the Q sr class. Thus we can write the equation for 
Q, 3 S k u in the form 

(2.16) PWS k U = Q j (S + 2) k Pu + Q U < 3j+ 9k-3 Q G Qsr- 

Suppose the function u solves the equation 

Pu = f. 

Commuting with all vector fields, we obtain equations for the functions u a . These 
can be written in the forms 

(2.17) Pu a = Qu <H + f< H =: F ai Q G Q sr , 

(2.18) au a = Qu< H + m 2 u< H +/< H =:G a , Q G Q sr , ReS z ad {r^)- 
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For F a and G a we have pointwise bounds of the form 

(2.19) \F a \ < J^(|tt 2 W< M- 6 | + |«<|a||) + ^(|V 2 «<|a|| + |V«< H |) + |/<|„||, 

(2.20) \G a \ < ^(|tt 2 u < H | + |«< w |) + JL(|vV w | + |Vu< w |) + |/< H |. 

As a general principle, we will use the latter equation to improve the bounds on 
u a away from r — (precisely for r > t), and the former near r = (precisely for 
r <C t). 

2.2. The weak local energy decay. The statement of the weak local energy 
decay property in Definition 11.31 includes the vector fields T but not S or ft. We 
remedy this in the following 

Lemma 2.6. Assume that the weak local local energy decay property (|1.9p holds. 
Then we also have 

(2.21) ll«<m|Usi < ||Vu<m+i(0)|| La + ||/< m +l|Ui5.. 



Proof. We use induction with respect to m. For m = the bounds (|1.10|) and 
(|2.21l) coincide. Consider now some m > 0, and a a multiindcx with |a| = m. If 
Z a contains only T derivatives then the bound for u a follows directly from ((1.1 OP . 
Else we factor 

Applying (fTTTOj) to n j S k u and using (|2.16[) we obtain 

\\Z a u\\ LEl < \\WS k u\\ LEi ,i 

< \\VQ?S k u(0)\\ H *+i + \\tt'S k f\\ LE i + 1 -* + \\Qu<3j+9k-3\\LE i + 1 -' 

< ||Vw< m+1 (0)|| L 2 + ||/< m+ i||i S . + ||(r)- 2 Vu< ro _i|U B . 

+ ||(r)" 3 M< m _i|| L£ ;» 

< ||Vu<„ l+ i(0)|| L 2 + ||/< m+ i||L£* + ||w<m-l||L£i 

which concludes our induction. □ 

2.3. The stationary local energy decay. Our first aim here is to include the 
vector fields S and in the stationary local energy decay bounds. A second aim is 
to derive a variation of the same bounds with different weights. 

Lemma 2.7. Assume that the stationary local energy decay property (|1.11[) holds. 
Then for all < to < t\ we also have 

(2.22) \\u< m \\ LE l[ toM ] < ^2 II Vu <™(^)IIl 2 + \\f<rn\\LE* + \\d t U< m \\ LE [ taM \, 

i=0,l 

respectively 

(2.23) ||V7i< m || L 2 < J2 \\{r)^u< m (ti)\\ L 2 + ||(r)/< m || L 2 + \\d t u< m \\ L 2. 

i=0,l 
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Proof. The proof of ()2.22|) is identical to the proof of Lemma \2. 61 and is omitted. 
We now prove (|2.23j) . We begin with the case m = 0, where we apply the classical 
method due to Morawctz. Assume first that we are in Case A. Multiplying the 
equation Pu = f by {xd x + l)u and integrating by parts we obtain 

Pu ■ {xd x + l)udxdt = ( ( \Vu\ 2 + 0((r)- 1 )\Vu\ 2 + 0((r)- 3 )\u\ 2 dxdt 



t JR 3 

0{{r))\Vu\ 2 + 0{{r)- 1 )\u\ 2 dx 

t=t 

Using Cauchy-Schwarz on the left and estimating the \u\ 2 terms by |Vu| 2 terms via 
Hardy inequalities we are left with 

LHSl^M(m = 0) < RHSgMim = 0) + || Vu\\ L * {r < R) + \\u\\ L 2 {r < R) 

for some fixed large R. Here the extraneous terms on the right are only measured 
for small r, as the large r contribution can be absorbed on the left. However, to 
bound them for small r we have at our disposal the bound (|2.22[) , whose right hand 
side is smaller than the right hand side of (|2.23[) . The same outcome is reached 
in Case B by inserting a cutoff function selecting the region {r 3> 1} in the above 
computation. 

To prove (|2.23|) we can use a simpler direct argument since there is no loss of 
derivatives on one hand and since we already have the bound (|2.22l) to use to 
estimate u <m inside a compact set. Precisely, for \a\ = m we have 

Pu a = f< m + 0((r)- 2 )Vu< m + 0{(r)- 3 )u< m . 

Then we apply the to — case of (|2 . 23[) for u a and sum over |a| < to. We obtain 

LHS^M Z RHS^M + l|Vu< m || L 2 (r < fl) + \\u< m \\mr<R) 

and the last terms on the right are estimated by (|2.22|) . □ 

3. The pointwise decay 

The strategy of the proof of our pointwise decay estimates is to itcratively im- 
prove the estimates via a two step approach. The two steps are as follows: 

(i) Use the properties of the fundamental solution for the constant coefficient 
dAlcmbcrtian □ via the equation (|2.18j) . This yields improved bounds for r ^> 1, 
but no improvement at all for r ~ 1. 

(ii) Use the stationary local energy decay estimates for the operator P in the 
region r<Ct. This allows us to obtain improved bounds for small r. The transi- 
tion from L 2 to pointwise bounds is done in a standard manner via Sobolev type 
estimates. 

3.1. The cone decomposition and Sobolev embeddings. For the forward 
cone C — {r < t} we consider a dyadic decomposition in time into sets 

C T = {T < £ < 2T, r < t}. 

For each CV we need a further double dyadic decomposition of it with respect to 
either the size of t — r or the size of r, depending on whether we are close or far 
from the cone, 



C T = [J Cfu (J C T 



■u 

KJi<T/4 l<U<T/4 
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where for R, U > 1 we set 

= C T n {R < r < 2R}, C% = C T n {U < t-r < 211} 

while for R = 1 and U = 1 we have 

Cf =1 = C T n {0 < r < 2}, C^ =1 = C T n {0 < i - r < 2} 

with the obvious change for C\ in Case B. By Cft and we denote enlargements 
of these sets in both space and time on their respective scales. We also define 

C< T ' 2 = |J Cf. 

R<T/4 

The sets and represent the setting in which we apply Sobolev embeddings, 
which allow us to obtain pointwise bounds from L 2 bounds. Precisely, we have 

Lemma 3.8. For all T > 1 and 1 < R,U < T/4 we have 
(3.24) 



respectively 



w \\L*{cKy 



(3.25) |ML-<c ?) < i E H^'HI L2( c ?) + S- E \\™Vw\\ LHe uy 

Proof. In exponential coordinates (s, p, w) with t = e s and r = e s+p the bound 
(|3.24l) is nothing but the usual Sobolev embedding applied uniformly in regions of 
size one. The same applies for p.25[) in exponential coordinates (s, p, uS) with t = e s 
and t-r = e s+p . □ 

Expressed in terms of the local energy norm, the estimate p.24|) yields 

Corollary 3.9. We have 

(3.26) IMUcc-/*)^ E W&rtwW^wy 

i+j<2 



3.2. The one dimensional reduction. A main method to obtain pointwise es- 
timates for u a is by using the positivity of the fundamental solution to the wave 
equation in 3 + 1 dimensions and the standard one dimensional reduction. 

For solutions to (|2.18[) with vanishing initial data, we may apply time translation 
invariance and assume without loss of generality that G a is supported in C = {r < 
t}. Then define 

2 

(3.27) H a (t,r) =Ell^G Q (t,r W )|! L2(s2) , 

o 

for G a as in (|2.18| . By the Sobolev embeddings on the sphere we know that 
\G a \ < H a . Let v a be the radial solution to 

(3.28) Uv a = H a , v a (0) = d t v a (0) = 0. 
Then we can compare 

\u a \ < v a . 
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We can rewrite the radial three dimensional equation (|3.28[) as a one dimensional 
problem 

{d 2 t - d 2 r ){rv a ) =rH a 

which has the solution 

1 f* f rJrt ~ s 
(rv a )(t,r) = - / / pH a (s,p)dpds. 

A Jo J\ r ~t+s\ 

Assuming that H a is supported in the forward cone t > r, this is rewritten as 
(3.29) r v a (t,r) = -f pH a (s,p)dsdp 

1 J Dtr 

where D tr is the rectangle 

D tr = {0<s-p<t-r, t-r < s + p <t + r}. 

In order to handle the contribution from the intial data in both cases and the fact 
that (|2.18p only holds outside of the cylinder R + x B(0, i?o) in Case B, we modify 
the above argument. The one dimensional reduction is only used to improve the 
bounds on u a for large r. Hence we can truncate the functions u a outside a large 
ball using a cutoff function Xout which is identically one for large r. Then we can 
redefine G a as 

With this choice for G a the bound (|2.20[) still holds. Moreover by truncating outside 
of a sufficiently large ball, Xoutu a has vanishing initial data. We can use the one 
dimensional reduction to obtain bounds for u a for large r, while for small r, we 
shall rely on the Sobolev-type estimates of Lemma 13.81 

3.3. An initial decay bound. Here we combine the above one dimensional re- 
duction with the local energy bounds in order to obtain an initial pointwise decay 
estimate for the functions u a . This has the form 

Lemma 3.10. The following estimate holds: 

(3-30) \u a \ < .y ~^ (||u<a+n|U^ + \\(r)f< a+ n\\LE*) 

(r) [t — r) 2 

Here and later in this section n represents a large constant which docs not depend 
on a but may increase from one subsection to the next. For the above lemma we 
can take n = 25 for instance, but later on it becomes tedious and not particularly 
illuminating to keep track of the exact value of n. We shall do so similarly for the 
enlargements of the cones Cj, . We shall not track, though we note that only a 
finite number will ever be required, each subsequent enlargement which is needed 
and shall allow the enlargement to change from line to line while maintaining the 
same notation. 

Proof. We assume that r> 1. For r < 1 we instead use directly the Sobolev type 
embedding (pT26]h 

For large r, we bound u a by v a and the function H ai using (|2.20|) , by 

\\(r) 2 H a \\ LE + \\(r) 2 SH a \\ LE < ||u< a +26||iBi + ||<r>/< a +26||iB.. 
Hence it remains to show that the solution v a to (|3.28j) satisfies 

(3-31) k| < ) ^~ r h W) 2 H a \\ LE + \\{r) 2 SH a \\ LE ). 

(r) [t — r) 2 
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We now prove (|3.31[) . The index a plays no role in it so we drop it. One can 
then estimate \rv\ as 

(3.32) \rv(t,r)\< [ pH(s,p)dsdp. 

JD tr 

We assume that r ~ t, as there is no further gain for smaller r in estimating the 
integral on the right. 

We partition the set D tr into a double dyadic manner as 

D tr = (J Dg, D* =D tr n{R<r< 2R} 

R<t 

and estimate the corresponding parts of the above integral. We consider two cases: 

(i) R < (t — r)/8. Then we need to use the information about Su. For any 
(s,p) G C, let 7s iP (t) be the integral curve corresponding to the vector field S, 
parametrized by time, satisfying 7 S;P (0) = (s,p). The fundamental theorem of 
calculus combined with the Cauchy-Schwarz inequality gives 

\H( ls , p m\ 2 <- f\H{ ls>p {r))\ 2 dT+- [ S \(SH)( 7stP (r))\ 2 dr. 
s Jo s Jo 

We apply this for (s, p) € D^ r and integrate. In the region D^ r we have p ~ R and 
\s — (t — r)| < R; therefore we obtain 

/ \H\ 2 dsdp<— [ \H\ 2 + \SH\ 2 dsdp 
Jd? t t-r J Br 

where 

B R = {{s,p) : R/8<p<8R}. 
Hence by Cauchy-Schwarz we conclude 

/ pHdsdp < —^(WHW^s^ + \\SH\\ L 2 (Br) ) 
Jd? t (t — r) 2 

< -\T{\\{r) 2 H\\ LE + \\{r) 2 SH\\ LE ). 
(t — r) 3 

The logarithmic factor in (|3.31|) arises in the dyadic R summation. We note that 
in the L 2 (Bji) norm above H is viewed as a two dimensional function, whereas the 
LE norm applies to H as a radial function in 3 + 1 dimensions. 

(ii) (t — r)/8 < R < t. Then we neglect SH and simply use Cauchy-Schwarz, 

/ pHdpds < R$(t - r) L 2\\H\\ L 2 (Bn) < R-^t-r^UrfHWLE. 

J Df T 

The dyadic R summation is again straightforward. □ 

3.4. Improved L 2 gradient bounds. The bounds obtained in the previous step 
for u a apply as well to Vu Q . However, u a and Vu a do not play symmetrical roles 
in the expressions for F a and G a . In particular, the weights that come with Vu Q 
are worse than the ones that come with u a . Hence, when we seek to reiterate and 
improve the initial pointwise bound (|3.30j) it pays to have better bounds for Vit Q . 
This is the aim of this step in the proof. Our dyadic L 2 gradient bound is contained 
in the next lemma: 
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Lemma 3.11. For 1 < U, R < T/4 we have 

(3.33) ||Vw||l2(c£) £ R^WMl l*(c$) + t ~ 1 \\ Sw \\lhc?) + r \\ Pw \\l^c^) 



respectively 

(3.34) l|vH| L2(C L f) < u^iWwh'iev) + \\ Sw h'(e$)) + T W Pw W l--{c^ 

Applied to u a this gives 
Corollary 3.12. For K [/, P < T/4 we have 

(3.35) l|Vu a ||z,2(C«) ~ R ~ 1 \\ u <a+n\\ L 2 { cR- ) + R\\f<a\\ L 2(cR) 

respectively 

(3.36) l|Vlia|| L 2 (c .c/) < ?/ -1 ||M<o+ re || L 2(c?) + T ll/<«lli=((5^) 
Applied to Vu Q we also obtain 

Corollary 3.13. Forl<^U,R< T/4 we have 

(3.37) l|V 2 u Q || L 2 (C H) < iJ -1 ||Vu< + re || L 2(cH) +- R IIW<a|lz,2 ( cR) 
respectively 

(3.38) ||V a «„|| La(0?) < C-^lVu^+nlli,^) + T||V/< a || L2(c5?) 

Proof of Lemma \3.11\ To keep the ideas clear we first prove the lemma with P 
replaced by □. We consider a cutoff function (3 supported in Cj$ which equals 1 on 
. Integration by parts gives 

(3.39) / f3(\S7 x w\ 2 - \d t w\ 2 )dxdt = [ Uw ■ (iwdxdt - \ [ (D(3)w 2 dxdt. 



To estimate Vw we use the pointwise inequality 

|VH 2 < M-^\Sw\ 2 + -^(\V x w\ 2 \dM 2 ) 
(t — ry t — r 

which is valid inside the cone C for a fixed large M. Hence 

/ [3\Vw\ 2 dxdt < [ . 1 R\Sw\ 2 + — —\up\w 2 + -^- [3\nw\\w\dxdt 
J J [t — r) £ t — r t — r 

where all weights have a fixed size in the support of [3. The function (3 can be 
further chosen so that |D/3| < r~ 2 . Then the conclusion of the lemma follows by 
applying Cauchy-Schwarz to the last term. The argument for C!f is similar, with 
the only difference that now we have |D/3| < £ _1 (£ — r)^ 1 . 

Now consider the above proof but with □ replaced by P. Then, given the form 
of P in (|1.2[) . p.3p . the relation p.39[) is modified as follows: 

/3{\ V x w\ 2 - \d t w\ 2 )dxdt = j Pw [Swdxdt - i J ((P + V)[3)w 2 dxdt 

+ J 0{r~ l )P\Vw\ 2 dxdt. 

The bound for (P + V)f3 is similar to the bound for Dj3, and one can easily see that 
the last error term is harmless. The proof of the lemma is concluded. □ 
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3.5. Improved L 2 bounds for small r. A very unsatisfactory feature of our 
first pointwise bound (|3.30[) is the r _1 factor which is quite bad for small r. Here 
we devise a mechanism which allows us to replace this factor by t . Our main 
bound is an L 2 local energy bound, derived using the stationary local energy decay 
assumption. We have 

Proposition 3.14. Assume that the problem (jl.ip satisfies stationary local energy 
decay bounds (jl.lip . Then the following estimates hold: 

(3.40) l|w<m|| Ls i (c <T/2 ) < r" 1 ||(r)u< m+ „|| Lsl(c <T/2 ) + \\f<m+n\\ LB „ (c <T/2 y 

<iT /2 

Proof. We first observe that we can harmlessly truncate u to C T . We will make 
this assumption throughout. Applying the stationary local energy decay estimate 
(|2.22j) for u< m we obtain 

\\u< m \\LEHc T ) ^ ||Vu< m (T)|| i 2 + \\Vu< m (2T)\\ L 2 

+ \\f<m\\LE'(C T ) + \\dlU< m \\LE(C T )- 

For the last term on the right, due to the support of u, we have 

\\dtU< m \\ L E(C T ) < 7p(\\ Su <m\\LE(C T ) + II rd r U< m \\ LE (Ct) ) 

< T-5(||( r )-i5 M < m || L2 + \\Vu< m \\ L ,) 

<T-*(||V5«< m ||£» + ||Vu< m ||i») 

where a Hardy inequality was used in the last step. Next we consider the time 
boundary terms. By the fundamental theorem of calculus and the Cauchy-Schwarz 
inequality, for any s > we have 

(3.41) |V^T,xO| 2 <^^(V U )( 7T ^(r))| 2 dr+^^ S |(5V U )( 7Ti:l .(r))| 2 dr. 

This holds uniformly with respect to s, so we have the freedom to choose < s < T 
favourably depending on x. Suppose we take s = T . Integrating the above estimate 
over \x\ < T for this choice of s and applying a change of coordinates yields 

/ \Vu< m {T)\ 2 dx <T~ l [ \Vu< m \ 2 + \SVu< m \ 2 dxdt. 

Jc T n{t=T} Jc T 

A similar bound holds for Vu< m (2T). Combining the last three estimates we obtain 

(3-42) WuKraWLE^iCx) ^ 1 1 f<m+n \\ LE* (C T ) + T ~ 1 1 1 ^ u <m+n \ \ L 2 (C T ) 

which brings us half of the way to the proof of p. 401) . 

We now consider the second expression on the right above, and we estimate it 
via the same argument as before, but using (|2.23|) instead of (|2.22l) . Indeed, (|2.23[) 
yields 

\\Vu< m \\ LH o T ) < ||<r)/<m||^(c T ) + \\(r)^u< m {T)\\ L 2 + \\{r)^-Vu< m {2T)\\ L 2 

+ \\d t U< m \\L 2 (C T )- 

The last term is controlled as above by 

||<9 t M< m || L 2 (CT) < r _1 (||5w< m || L 2 (c . ;r) + ||rV-u< m || L 2 (C;r) ) 

ST~* (\\(r)Su< m \\ LE i(c T ) + \\{r)u< m \\ LE i(c T ))- 
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The initial and final terms are estimated by a more careful use of (|3.41j) . Namely, 
we integrate p.4ip with s(x) = (r) 3 T*. This yields 



J (r)\Vu< m (T)\ 2 dx < T-* £ (r)i\Vu< m \ 2 +T 



(r)i\SVu< m \ 2 dx dt 



which implies that 

\\(r)^u<m(T)\\ L 2 <T-i||(r)iVu< m ||L»(c r )+r-*|Kr)5u< m || LB i C c T ). 

Combining the last three bounds leads to 

||Vu< ln || L 2 (cv) < ||(r)/< m+ „|| L 2 (CT) + T~i\\(r)u< m+n \\ LE i {CT) 

+ T-3||(r)3Vu< m || L 2 (CT) . 

We can discard the last term on the right by absorbing it into the left hand side 
for r <C T and into the second right hand side term for r ~ T. Thus we obtain 

(3.43) ||Vw< m || L 2 (CT) < \\(r)f< m+n \\ L 2 {CT) + T~ || (r)u< m+n \\ LE i {CT y 

Finally, the bound (j3~4T)| is obtained by combining ((3~4^|) and (fg~S5|) . □ 

3.6. Improved pointwise bounds for u a for small r. Here we use Sobolev 
type inequalities to make the transition from L 2 local energy bounds to pointwise 
bounds for small r. Our main estimate has the form 

Proposition 3.15. We have 
(3.44) 

3 

||u< m || LOO(c < r/3) + ||(r)Vu< m || ioc(c <T/2 ) < T-*||u< m+n || L£ , (< 5<T/ a) 

+ T ~~- (ll/<m+n|| LE *(C< T/ ' 2 ) + V /<m+n II LE(C< T ' 2 ) ) ' 

Proof. In view of the Sobolev bound (|336|) in Corollary [3JB the estimate (|3~H|) 
would follow from 

W u <™\\lehc< t/2 ) + \\^ Wu < m hEHc< T/2 )^ T ' 1 ^ u < m + n h E (c< T/2 ) 
(3.45) 2 

+ ll/<™+n|l L£ ;* ( c< T / 2 ) + IK r ) V f<™+n\\ LE( c<m-y 

To prove this we begin with the local energy bound (|3.40[) and add to it corre- 
sponding bounds for the second derivatives of u< m . For that we use f|3 .3T|) : applied 
uniformly in dyadic regions Cft with R < T it yields 

||{r)V 2 u< m || jL£7(c ,<T-/2 ) < \\u< m+n \\ LE1(c < T/ 2 ) + \\{r) 2 V .f< m+n \\ LE[c <T/2 y 

Combined with (|3.40p . this leads to 

LHS®A5$ < T~ 1 \\(r)u< m+n \\ LE1[d< T / 2 ) + \\.f< m +n\\ le , {c <t/2 ) 

(3.46) 

+ Vf<m+n\\ LE{c <T /2y 

It remains to make the transition from (|3.46[) to (|3.45[) . This is done via (|3.35[) 
applied uniformly in dyadic regions Cft with R<T, which yields 

\\( r ) u <m\\ LE i ( c< T / 2 ) ^ \\ U <™+™\\lE(C< t/2 ) + IK r ) 2 ^<™+«llLB(C< T/2 )- 

Hence (|3~45)) follows. □ 
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3.7. Improved pointwise estimates for Vit Q near the cone. Here we convert 
the improvement in the L 2 bounds for Vu< m given by Lemma 13.111 near the cone 
into a similar L°° bound. This shows that the pointwise bounds for Vit< m are 
better than those for u< m by a factor of (t — r) _1 . 

Proposition 3.16. We have 

U\\Vu< m \\ L °e(C") ~ ll u <m+n|li-(C^) + J1 " 5C/5 ll/<™+™llL2((5i!) 

(3.47) 

+ T-5C/i||V/< m+ „|| i2(d?) . 



Proof. The proof is similar to the proof of Proposition 13. 15l In view of the Sobolcv 
inequality (|3.25[) . the bound (|3.47[) would follow from 

II Vu< m || + C/ 2 ||V 2 U< m || L 2 (c u ) < \\u< m+n \\ L 2 (c u ) +UT\\f< m+n \\ L 2 (c u ) 

+ [/ 2 r||v/< m+n || L 2 (c .c/). 

To prove the bound on Vu< m we use (|3.36[) . which shows that 

U\\ Vw< m || L 2( C C() < ||w< m +, l || i 2( ( f.C/) + UT\\f< m+n \\ L2 (guy 

To prove the bound on V 2 it< m we use (|3.38[) to obtain 

/7 2 ||V 2 M< m || i 2( C c/ ) < t/||Vu< m+ „|| L2(( =,u ) + U 2 T\\\7f< m+n \\ L2 ^uy 

□ 

3.8. Gradient bounds associated to the first decay bound. Here we start 
with the bound (|3.30p for u< m and improve it for small r, as well as complement 
it with gradient bounds. 

Lemma 3.17. The following pointwise estimates hold: 
(3-48) |w<m|^Ci— |VM< m |<Gi- 



t(t-r)i (r)(t-r)i 
where 



1 „1 1., „ m 1 _ 3 ,. T 3 , 



Cl =||u< m + „|i Li Jl+SUpT2i?2[/2||/< m+n || i2(c , R|t/) +T 2i?2[/2||V/< m+ „|| i2(c ,H, C / ) . 

-R, U 

Here for brevity C^' U stands for either Cft or , with the natural conven- 
tion that R ~ T in and U ~ T in . The proof is a direct application of 
Propositions 13.151 13.161 using p.30[) as a starting point. 

3.9. The second decay bound. 

Lemma 3.18. The following decay estimate holds: 

(3.49) | u < m | < C 2 i#^, |V U < ro | < C 2 l06{t - r) 



t(t-r) ' 1 '(r)(t 



r 



2 



whe 



C 2 = \\u< m+n \\ LE i +SUpTi?2j72 ||/<m+n|| z ,2( C R,c- ) + -R 2 E/ 2 || V/< m+n || jj^b.Bj . 
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Proof. By the Sobolev embeddings of Lemma 13.81 we have 



^ ~ t(r)*\t-r) C2 - 



Also by (|3.48p we have 



F<m+el ^ — |Vu< m+6 | ^ — — -j^- 

t{t — r) 2 (r)(t — r) 2 

Hence for the functions G< m we obtain 

|G ^ m| ~ (r)3(Lr) C2 - 

Computing via the one dimensional reduction, this leads to a bound for u< m of the 
form 

\U<m\ < C 27T7 7 T 

(r)(t-r) 

which is comparable to (|3.49p near the cone, but it is weaker for r < f. Then 
the small r bound for u< m and the bound for Vu< m arc obtained from Proposi- 
tions [37TH1 GEES (we need to increase n appropriately at this stage). □ 

3.10. The third decay bound. 

Lemma 3.19. The following decay estimate holds: 

(3.50) c s i|i^>, iv^isftj*^ 

where 

C.3 = ||li<m+n|UBi + SUpT 2 i?5[/5||/< m+Jl || L2(c ^ iLf) + TR^lji || V/< TO+ra || i2j - c H,l7). 



Proof. As before, the main step in the proof is to obtain a pointwise bound for it< m 
which coincides with (|3.50[) near the cone, 

(3-51) |«< ro | < ^. { *~^ C 3 . 

Once this is done, the full estimate ()3.50j) follows easily by a direct application of 
Propositions I3.16[ 13.151 However, at this stage there is a new twist, namely that 
the one dimensional reduction no longer suffices for the proof of (|3.5ip . 
The pointwise bound for f< m has the form 

l/<™l £ t 2 {r) 2 {t _ r) C *- 

Also by (|3.49[) we have 

i | < log(t - r) log(t - r) 

\u< m+ s\ < -^37)-^, |V W <„ 1+6 | < (r)(t _ r)a g 3 - 

We use the wave equation for u< m given by (|2.18p . and rewrite G< m in the form 
G< m = f<m + au< m+6 + d t (bu <m+6 ), a e S z (r~ 3 ), b G S z (r~ 2 ). 
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Here we can confine ourselves to dt derivatives in the last term because for any S 
and ft component we gain a factor of and include it in the second term. Wc 
split G< m into two parts, 

G<m = G< m + C<m 

with 

G <m = f<m + au< m+6 + d t ((l - x(t, r))bu <m+6 ), G< ln = <9t(x(t, r)bu <m+6 ) 

Here x is a smooth cutoff selecting the region t — r < i. 

The function G 1 <m contains the part of G< m which has good pointwisc bounds, 

, r i | < log(*-Q 

Computing via the one dimensional reduction, this gives the pointwise bound (|3.51[) 
for the corresponding part u\ m of u< m . 

Next we prove the same bound for the output u< m of G 2 <m . Taking absolute 
values and applying the one dimensional reduction does not work, as it misses a 
cancellation due to the presence of the derivatives. Instead we do a more precise 
computation 

Lemma 3.20. Consider a smooth junction f supported in {| < r < t} so that 
(3.52) \f\ + \Sf\ + \Qf\ + (t-r)\d r f\ < 



t 3 (t - r) \og 2 (t - r) 



Then the forward solution u to 

satisfies the bound 
(3.53) 



□u = dtf 



< 1 



t(t- r y 



We note that if (|3.52p is replaced by 
(3.54) |/| + |V/| + \Sf\ + \nf\ + {t- r)\d r f\ < 1( 



t 3 (t-r) 

then (|3.53| is trivially replaced by 
(3.55) H< l0g3( *- r) 



t{t-\ 



) 2 

which suffices to conclude the proof of Lemma 13.191 



Proof. The function u is expressed in the form u = d t v with v the forward solution 
to Ov = f. Via the one dimensional reduction applied to v. Vw Civ, Sv and 
(tdi + Xidt)v we obtain 

M + |V«| + \Sv\ + \flv\ + J^Ktdi + Xi dt)v\ < 

tit — r) 

i 

where the main contribution comes from the cone. The above left hand side domi- 
nates (t — r)dtv; therefore the proof of the lemma is complete. □ 

□ 
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3.11. The fourth (and final) decay bound. Here we reiterate one last time to 
remove the logarithms in (|3.50[) and establish the final bound, which concludes the 
proof of the Theorem. 

Lemma 3.21. The following decay estimate holds: 
(3.56) |u< m | < C V /4 1 , 2 , |Vu< m | < C 4 



t(t-r) 2 ' 1 - ' ~ *<r)(i-r) 3 

where 

Ci = \\u< m ^\\ LE i+supY,T 2 RiU^\f< m+n \\ L2{C R,u ) +TRiui\\Vf<r, 

T R,U 

Proof. The argument is similar to the previous step, but with some extra care in 
order to avoid the logarithmic losses. The main goal is again to obtain a pointwise 
bound for u< m which coincides with (|3.50[) near the cone, 

(3-57) |u< m | < TT7 - — -2-C4 

(r)(t-r) 2 

after which the full estimate (|3.50j) follows from Propositions 13.161 13.151 
The pointwise bound for f< m still has the form 

l/<™l ~ t 2 {r) 2 {t _ r) C ^ 

but now in addition we have dyadic summability with respect to R and U. Also by 
(|3.50p we have 

i log 3 (t-r) ^ | < log 3 (t-r) ^ 

We split G< m = G 1 <m + G\ m exactly as before. 

For G 1 <m we use the one dimensional reduction, based on the bound 

\G\ra\ ^ \ f<m\ + J^TZT^3 C ^ 

which gives the pointwise bound p.57p for w< m . The dyadic summability for / 
causes the absence of logarithms in the bound for the contribution of /. The 
contribution of the second term is of the order of 

log 3 (^-r) 
(r){t-r)* 

where we have a full extra power of (t — r) available to control the logarithms. 
Finally, the contribution of G 2 <m is controlled by Lemma [3.201 □ 

4. Perturbations of Kerr spacetimes 



We now present an application of Theorem II .51 to general relativity. We are able 
to recover Price's Law not only for Schwarzschild and Kerr spacetimes, but also for 
certain small, timc-dependent perturbations thereof. We begin by presenting the 
results obtained in [45] and (66] for the Schwarzschild and Kerr metrics. We continue 
in 14.21 with a proof of stationary local energy decay estimates for perturbations 
of Schwarzschild; while the perturbations are required to be small, no decay to 
Schwarzschild is assumed. This result applies as well to small perturbations of 
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Kerr with small angular momentum. Finally, in 14.31 we establish weak local energy 
decay estimates for small perturbations of Kerr; here, we essentially require a 
decay rate of the perturbed metric to Kerr. 

4.1. The Schwarzschild and Kerr metrics. The Kerr metric in Boyer-Lindquist 
coordinates is given by 



ds = gttdt + gt(f>dtd(f) + g rr dr 

where t G M, r > 0, (<f>,0) are the spherical coordinates on § 2 and 

A-a 2 sin 2 6> 2Mrsin 2 6» p 2 

gtt = o ' 9t<j> = -2a - , g rr — — 

p l p 2 - A 

(r 2 +a 2 ) 2 - a 2 Asin 2 6> . , . , 

gu = t, sin 9, gee = p 

P 2 

with 

A = r 2 -2Mr + a 2 , p 2 = r 2 + a 2 cos 2 9. 

Here M represents the mass of the black hole, and aM its angular momentum. 
A straightforward computation gives us the inverse of the metric: 

„ (r 2 +a 2 ) 2 -a 2 Asin 2 u 2Mr rr A 

9 tt = -- 9~r , g^ = -a 



p 2 A ' » p 



2 ' 



4>4> = A - a 2 sin 2 9 ee 1_ 

9 p 2 Asm 2 9 ' 9 p 2 ' 

The case a = corresponds to the Schwarzschild space-time. One can view 
M as a scaling parameter, and a scales in the same way as M. Thus M/a is a 
dimensionless parameter. We shall subsequently assume that a is small a/M <C 1, 
so that the Kerr metric is a small perturbation of the Schwarzschild metric. One 
could also set M = 1 by scaling, but we prefer to keep M in our formulas. We let 
.9s i 5k denote the Schwarzschild, respectively Kerr metric, and ds, Qk denote the 
associated d'Alembertians. 

In the above coordinates the Kerr metric has singularities at r — on the equator 
9 = tt/2 and at the roots of A, namely r± = M ± \J M 2 — a 2 . The singularity at 
r = r + is just a coordinate singularity, and corresponds to the event horizon. The 
singularity at r = r_ is also a coordinate singularity; for a further discussion of 
its nature, which is not relevant for our results, we refer the reader to 130] . To 
remove the singularities at r = r± we introduce functions r* , v + and </>+ so that 
(see [30]) 

dr* = (r 2 + a 2 )A~ 1 dr, dv + = dt + dr*, d<f> + = d(f> + aA' 1 dr. 

We call v+ the advanced time coordinate. The metric then takes the Eddington- 
Finkelstein form 

2Mr 

ds 2 = - (1 Tr-)dv% + 2drdv + - 4ap~ 2 Mr sin 2 9dv + d(j) + - 2a sin 2 9drd<j) + 

P 

+ p 2 d9 2 + p- 2 [(r 2 + a 2 ) 2 - Aa 2 sin 2 9} sin 2 9d<j>\ 

which is smooth and nondegenerate across the event horizon up to but not including 
r = 0. 

In order to talk about perturbations of Kerr we need to settle on a suitable co- 
ordinate frame. The Boyer-Lindquist coordinates are convenient at spatial infinity 
but not near the event horizon while the Eddington-Finkelstein coordinates are 
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convenient near the event horizon but not at spatial infinity. To combine the two 
we replace the (i, <f) coordinates with (t, cf>) as follows. 
As in [35] and [56], we define 

t = v+- n(r) 

where p is a smooth function of r. In the (t,r,cj) + ,0) coordinates the metric has 
the form 

ds * = (1 _ ^L) d P + 2 (l - (1 - ^fV(0) dtdr 

/ 2Mr \ 

— 4ap- 2 Mr sin 2 9did<f>+ + \2p\r) - (1 - — — ){p\r)) 2 j dr 2 

- 2a0(l + 2p- 2 Mrfi'(r)) sin 2 drd<j)+ + p 2 d0 2 
+ p~ 2 [(r 2 + a 2 ) 2 - Aa 2 sin 2 9] sin 2 0d<p 2 + . 

On the function p we impose the following two conditions: 

(i) p(r) > r* for r > 2M, with equality for r > 5M/2. 

(ii) The surfaces t = const are space-like, i.e. 

2 Mr 

U?(r)>0, 2 - (1 - — )/x'(r) > 0. 

For convenience we also introduce 

= C(r)0+ + (1 - CW)0 

where C is a cutoff function supported near the event horizon and work in the 
(t, r, <p, 9) coordinates which are identical to (t, r, <fr, 9) outside of a small neighbor- 
hood of the event horizon. 

Given r_ < r e < r+ we consider the Kerr metric and the corresponding wave 
equation 

(4.1) D K u = / 
in the cylindrical region 

(4.2) M = {t > 0, r > r e } 
with initial data on the space-like surface 

(4.3) £- = Mn{t = 0}. 
The lateral boundary of M , 

(4.4) £+ = M n {r = r e }, = S+ n {t < t < h} 

is also space-like and can be thought of as the exit surface for all waves which cross 
the event horizon. This places us in the context of Case B in section 11.11 The 
choice of r e is not important; for convenience one may simply use r e — M for all 
Kerr metrics with a/M <C 1 and small perturbations thereof. 

A main difficulty in proving local energy decay in Schwarzschild/Kerr space- 
times is due to the presence of trapped rays (null geodesies). In the Schwarzschild 
case, this occurs on the photon sphere {r = 3A/}. Consequently the local energy 
bounds have a loss at r = 3M. To localize there we use a smooth cutoff function 
Xps( r ) which is supported in a small neighborhood of 3M and which equals 1 near 
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3M. By Xps(r) we denote a smooth cutoff that equals 1 on the support of Xps- 
Then we define suitable modifications of the LE 1 , respectively LE* norms by 

3M 

IMIlbi = \\d r u\\ LE + ||(1 — )Vu\\le + \\r~ u\\le 

\\I\\lE* = 11(1 - Xps)f\\LE* + ||Xps/!l_ffi + (l-Mi)L2 

With these notations, the local energy decay estimates established in [45] have the 
form 

Theorem 4.1. Let u solve Dsti = / in A4. Then 

(4-5) \\u\\ LE i +sup||V U (i)|| L 2 < \\Wu(0)\\ L , + \\f\\ LES . 

As written there is a loss of one derivative at r — 3M. This can be improved to 
an e loss, or even to a logarithmic loss, see [IS], but that is not so relevant for our 
purpose here. 

In the case of Kerr, the trapped rays are no longer localized on a sphere. However, 
if a/M <C 1 then they are close to the sphere {r = 3M}. 

We now recall the setup and results from [66] for the Kerr spacetime. Let r, £, $ 
and 9 be the Fourier variables corresponding to t, r, <j) and 0, and 

p K (r, 9, t, C, $, 6) = g tt T 2 + 2 5 *^r$ + g^® 2 + g rr f + g ee G 2 

= g u (T - n(r, 9, $, 0))(t - r 2 (r, 9, £, $, 9)) 

be the principal symbol of Dk- Here n, r 2 are real distinct smooth 1-homogcneous 
symbols. It is known that all trapped null geodesies in r > r + satisfy 

(4.6) R a (r,T,$)=0 

where 

R a (r, t, $) = (r 2 + a 2 ) (r 3 - 3Mr 2 + a 2 r + a 2 M)r 2 - 2aM(r 2 - a 2 )r$ - a 2 (r - M)$ 2 . 

Let r a (r, $) be root of (|4.6jl near r = 3M, which can be shown to exist and be 
unique for small a. Then define the symbols 

c i (r,e,e,*,e)=r-r (T i ,*), i = l,2 

and the associated space-time norms: 

||u||i aot = ||c i (D, a ;)u||i a + |H|| r - 1 

The replacements of the LE and LE* norms are 

IMIlbi = \\Xps(Dt-T 2 (D,x))c 1 (D,x)xp S u\\ 2 L 2 



+ \\x P s(D t - t 1 (D,x))c 2 (D,x)x P su\\ 2 L 2 

,2\ 



+ \\d r u\\ LE + ||(1 - x 2 )^u\\ LE + Wt^uWle 

ll/IU-EJc =11(1 - Xps)f\\LE* + \\Xpsf\\ Cl L^nc 2 L 2 - 

Then the main result in |66j asserts that 

Theorem 4.2. Let u solve O k u = f in M.. Then 

(4-7) \\u\\ LE i + sup \\Vu(t)\\ L 2 < ||V«(0)|| La + \\f\\ L E k . 

K v>0 
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This is the direct counterpart of (|4.5[) . which corresponds to t\ = — T2 and 
ci = C2 = r — 3M. Again, the loss of one derivative can be improved to an e loss, 
or even to a logarithmic loss, see |67j . 

4.2. Stationary local energy decay for small perturbations of Schwarzs- 

child. Here we consider a Lorcntzian metric g in A4 which is a small perturbation 
of Schwarzschild expressed in the (i, r, <j>, 0) coordinates. Our main result is as 
follows: 

Theorem 4.3. Let g be a Lorentzian metric on M and u a smooth function in M. 

a) Let Xps be a smooth cutoff function which selects a small neighbourhood of the 
photon sphere {r = 3M}. If 

(4.8) \d a [g^-(g K )^}\ <e(r)-^-, < |a| < 1 

for a small enough e, then the stationary local energy bound holds for u: 

(4-9) \\u\\ LE1[ioA] < ||Vu(to)|| L 2+||Vu(?i)|| L 2+||a sU || iJ 5, [t ~ 0ifl] + ||x ps 5 t -u|| iJ5[t - 0it - l] . 

b) If in addition 

(4.10) I^VI <r- x ~, l<|a|<fc+l 
then (jl.lip also holds. 

Proof, a) The argument is based on the computation in |45j . which only requires 
the use of vector fields. We will first prove the estimate in the Schwarzschild case, 
but do it in such a way so that the transition to g is perturbative. 
Let X be a differential operator 

(4.11) X = b{r)d r + c{r)d~ t + q{r) 

for some smooth functions b,c,q : [r e , oo) — > R with c constant outside a compact 
region and b, q satisfying 

, x \d?b\ < cr-W 

4 - 12 1 , 

\d?q\<c a r~ 1 -\ a \. 

Let -M[t tji = {to < t < ii, r > r e }, and let dVs = r 2 drdtduj denote the 
Schwarzschild induced measure. It was shown in [45] that one can find X as above 
so that 

(4.13) / Q s (t)dt = - I D s u- XudVs- BDi? s [u]|~^- BDR s [u]\ r=r 

Jin J M. r r_ jr. i 



with 



tl 

to 



" to 

where 

HIl^. = \\(r- mxVu\\l 2 + ||r- 2 a rU ||| 2 + Hr-^Hi, + ||(1 - X 2 )V^||i 2 
and the boundary terms satisfy 

BDR s [u]\ i=ii « ||V«(t 4 )lli=i. * = 1 > 2 

1 [to.'il' 
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Comparing the LE^ w norm we obtain from this computation with the LE 1 norm 
which we need, one sees that two improvements are necessary, one near r = 3AI 
and another for large r. 

The improvement for large r is a consequence of the fact that for large r one can 
view the Schwarzschild metric as a small perturbation of the Minkowski metric. 
Precisely, from the estimate [49] (2.3) (see also [48]) we have for large R 

(4-14) \\x>Ru\\LE 1 [t ;t 1 \ % WAleI^M + \\x>R a su\\ LE[io:il] + ||Vu(f )||i2. 

The similar bound for the metric g is also valid. 

To gain the improvement for r close to 3M we add a Lagrangian correction term. 
Precisely, a direct integration by parts yields 

a s u- X 2 ps udV s = I x 2 ps9s" d n ud v u +( a sX 2 p s)u 2 dV s + I xl s 9s^ t udx 

r„ r. , JMti. z. , J Mz 



t=t 

Since dt is timclikc in the support of Xps, we can write the pointwise bound 

lx ps Vu| 2 < x 2 ps 9s 'd mu ud v u + C\x ps d t u\ 2 

for some large constant C . Then the previous identity yields 
(4.15) 

Wx P sVu\\ 2 L2(M[ia j < f a s u-x 2 ps udV s +C\\x P sdtu\\ 2 L2{M[io h]) + \\Vu(U)\\ L 2. 

Combining (|4~T3)) . (|47l4|) and (|4T5)) we obtain 
(4.16) 

lkll!sHM r{ x ,) + Mhi(&- ) + < - f a sU - X lU dV s 

+ l|Vw(t )Hi2 + C \\Xpsd t u\\ 2 L 2 {M[igh]) + \\x>RO s u\\ 2 LEHM[iah]) 

where 

X, = X + S X 2 ps (r) 

with a fixed small constant S. At this point, the desired conclusion (|4.9j) for the 
Schwarzschild metric follows if we estimate the integral term by ||C1su||l.e*|I u II.l.e 1 - 
It remains to show that a similar estimate holds with □§ replaced by O g , This 
substitution is easily made in the last term on the right by performing a similar 
substitution in (|4.14p . Consider now the difference in the integral term, 



D= f (n s -n g )u-X lU dV s . 

J M<z~ i-.i 



To estimate this we use the bound (|4. 8[) to write 

□ s - U g = d^gg - g^)d v + 0{er- l -)V. 

Then we integrate by parts in a standard manner. Using also Hardy type inequal- 
ities we obtain 

\D\ <e(J2 HV«(t,-)ll£* + ||«||^ + ) + \\u\\ 2 LEl{M j). 

Hence (|4.16[) for D g follows, and the proof of the stationary local energy bound 
(|4.9[) is concluded. 



26 



JASON METCALFE, DANIEL TATARU, AND MIHAI TOHANEANU 



b) The proof follows closely that of Theorem 4.5 in [66], but for the sake of 
completeness we include it here. The result will follow by induction on k. The case 
k = is part (a) of the theorem. We will prove the case k = 1, and the rest follows 
in a similar manner. 

We need to estimate || V 2 u\\ LE ^ o ji, We already control || V%it|| iB rg ji, there- 
fore it remains to estimate the second order spatial derivatives. We write 

O g u = Lid^u + L 2 u 

where L\ is a first order operator and L 2 is a purely spatial second order operator. 
Then we have 

\\ L 2u\\ LE[ t a ; tl] < ||n s w||LB[t ,ti] + \\ L l d t u \\LE[t M\ 

which is favorable where L 2 is elliptic. But L 2 is elliptic wherever is time-like. 
Since g is a small perturbation of g$, this happens everywhere outside a small 
neighbourhood of r = 2M . Thus we have the elliptic estimate 

(4.17) ||.wV 2 M || LB[foii] < \\U g u\\ LE[ioM + \\Vd iU \\ LE[ioM + \\Vu\\ LE[ioM 

where Xout selects the region {r > 2M + 5}. 

It remains to estimate || ^ 2 Xeh u \\h 2 where Xeh is a smooth cutoff function which 
selects the region {r < 2M + 26} near the event horizon. The function v = Xeh u 
solves the equation 

U g V = hl= Xehf + Pg, Xeh] 

where the second term on the right is controlled in H 1 via (|4.17|) . Hence the 
conclusion of part (b) of the Proposition would follow from the following 

Lemma 4.4. Let M. e h = {r e < r < 2M + 3S, t > 0} with a fixed small 5. Let g be 
an 0(e) perturbation of gs in C m+1 (Ai e h) with e sufficiently small. Then for all 
functions v with support in {r < 2M + 3(5} we have 

(4-18) \\Vv\\ Hm{Meh ) < ||V«(0)|| H «. + \P g v\\ Hm(Meh ) 

The similar estimate holds in any interval [to,ti\. 

Proof. This is an estimate which is localized near the event horizon, and we will 
prove it taking advantage of the red shift effect. In microlocal terms, the red shift 
effect is equivalent to exponential energy decay along the light rays which stay on 
the event horizon, and small perturbations thereof. But for this estimate, these are 
all light rays of interest. All others exit the domain M. e h in a finite time. 

We begin with a simplification. If e is small enough then for r < 2M — 5 the r 
spheres are uniformly time- like, therefore we can use standard local energy estimates 
for the wave equation to reduce the problem to the case when r e = 2M — 25. 

For m = the above bound follows from part (a) of the Proposition. For m = 1 
we commute n g with the vector fields c^, H and d r . We have 

[d h a g ]=0(e)Q 2 , [n,a g ]=0(e)Q 2 

for some second order partial differential operator Q 2 with bounded coefficients. 
Hence applying (|4.18|) with m = to d^v and fif we obtain 

(4.19) \\&o\\m(Meh) + \\ d A\H^{M eh ) ^ \\Hm{M ch ) + ^Mhhm^)- 
We still need to bound d r v. For that we compute the commutator 

(4.20) p g ,d r ) = ~(d r g r s r )d 2 r + 0(e)Q 2 + N 2 
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where N2 stands for a second order operator with no d% terms. The key observation 
here is that 7 = <9 r <?g r > near r = 2M. We can now write 

(□ B - 7iX)a r « = d r h + (0(e)Q 2 + N 2 )v, 7l > 

with N 2 as above and most importantly, a positive coefficient 71. We recall here 
that X looks like — d r near the event horizon. Because of this the operator 

B = U g - lx X 



satisfies the same estimate (|4.9j) as n g for functions supported near the event hori- 
zon. To see this it suffices to examine (|4.16j) with Dg replaced by D g . Since X = X\ 
near the event horizon, it follows that the contribution of 71 X has the right sign 
and can be discarded. Hence we obtain 

^ 21 , \\d r v\\m(M eh ) < \\h\\m(M ch ) + 4 v \\H*(M eh ) 

+ \\tov\\Hi(M eh ) + \\dtl>\\Hi(M eh ) + \\v\\Hi(M sh ) 

where the last three terms account for the effect of N 2 . Combining the bounds 
(|4T9| and (|4~2"Tj) we obtain 

IMIff^A^) ^ WHH^Mih) + 4 V \\H*(Meh) + \\v\\Hl(Meh)- 

If e is sufficiently small then the conclusion (|4.18l) follows for m — 1. The argument 
for m > 1 is similar. □ 

□ 

4.3. Local energy decay for small perturbations of Kerr. Here we consider 
small perturbations of a Kerr space-time with small angular momentum, \a\ <C M. 
Our main result is as follows: 

Theorem 4.5. Let g be a Lorentzian metric on M. and u a function in M. solving 
Z g u = f. 

a) Assume that g satisfies (|4.8[) and decays to c/k near the photon sphere, 

(4.22) X P s\d a [g^-(g K )^}\ <c a (t), < |a| < 1 

where c a £ L~ (in particular, we can take c a = (t)^ 1 ^ ). Then the weak local energy 
estimate holds: 

(4.23) \Wf H1(J1+) + sup \\\7u(t)\\h + \\u\\l E , < \\Vu(0)\\h + \\f\\ 2 LE , 

R i>0 

b) Assume in addition that (|4.10[) holds. Then 

(4.24) ||u||* + . + sup ||Vu(t)||^ + \\u\\ 2 LEl , k < ||Vu(0)||^ + 

Rl ( >0 K K 



mi 



d thus (ll~9l) holds 



Proof. On any fixed compact time interval we have uniform energy estimates. Elim- 
inating a compact time interval, we can assume without any restriction in generality 
that the intcgrability condition on c(t) is strengthened to 



(4.25) / c(t)dt<e, \c(t)\<e 







where (|4.8j) was also used. 

a) The proof of (|4.23|) is similar to the proof of Theorem 14.31 but it requires 
the use of pseudodiffcrcntial operators. Let us start by recalling the idea behind 
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the proof of Theorem 14.21 It is shown in [66] that there exists a pseudodiffcrential 
operator Si of order 1 that satisfies the following: 

a) Si is a differential operator in t of order 1 . 

b) Si = X + XpsS w Xps, where X is defined as in (|4.11[) and s £ S 1 . 

c) Let .Mmjgi = {0 < i < to, r > r e } and dV^ = p 2 drdtdui denote the Kerr 
induced measure. Then one has 

(4.26) 

Q K (t)di = - f (a K u)(S lU )dV K - BDR K [u]\lZ - BDR K [u]\ r=r 



with 



Q^(t)dt>\\u\\i EkJoM 
\\ u \\' 2 LE^ m = \\Xps{Dt~T 2 {D 1 x))c l {D 1 x)x P sAv2 

+ \\Xps( D i -n(D,x))c 2 (D,x)x ps u\\ 2 L 2 

+ \\r- 2 d r u\\l 2 + llr-^Hl, + ||(1 - X 2 ps )Vu\\h 

and the boundary terms satisfying 

BDR K [u]\ i=h « ||V M (t i )|| L2 

(427) BDRK[u]\ r _ r *\\u\\ 2 + 

\r-r e a (^ [0tto] ) 



Note that conditions a) and b) guarantee that the boundary terms are well-defined 
after integrating by parts. 

The same reasoning as in Theorem 14.31 leads to the counterpart of (|4.16p . namely 



Ml 2 



lW(Mo, t - 0) ) + ll Vu ll^(^ fo p + ll v ^)lli 2 < 

(4.28) f °_ 

- / D K u ■ Sm dV-K + \\Vu(t )\\ L 2 + \\x>R a Ku\\LE*(M l0 j o] )- 

Here we seek to replace the Kerr metric by g. As discussed in Theorem 14.31 the 
bound (|4.14p holds as well for the metric g, therefore the last term is not an issue. 
Hence it remains to consider the difference 



D = (D K - O g )u ■ Siu dV K 

JM,n ST., 



We split Si = Xi + S ps where 

X\ = X — XpsXxps, S ps = XpsXxps + XpsS w Xps- 

Thus X\ is a first order differential operator which is supported away from the 
photon sphere. Correspondingly we split D = D\ + D ps . For D\, integration by 
parts using (|4.8[) leads to 



\Di\ < e(||Vu(0)||i a + ||Vu(io)||i 2 + ||«||^ 1(E+ j + M 2 leUm, f ,))• 

( [o,t ]' K [°i*o]' 

Here it is essential that the outcome of the integration by parts is supported away 
from the photon sphere {r = 3M}, where the LE^ and LE 1 norms are equivalent. 
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To estimate D ps we need to use the stronger bound (|4.22|) . Then near the photon 
sphere we can write 

1 



□ K - □„ = -=^VW« - <T)dv + 0(c(t))V 
V\9k\ 

where we also have gg - g^ v = 0(c(i)) and V(s&" - g^ v ) = 0(c(t)). Thus inte- 
grating by parts we obtain 

\D ps \ < f x ps cW(|V U | 2 + | M | 2 )^K + c(0)||V M (0)||i 2 + C (i )||V U (t )||| 2 . 

jM lo,t ] 

By ({4725]) it follows that 

\D ps \<e( sup \\Vu(t)\\l* + \H\IeUm rot - 

Applying the bounds for D\ and D ps , we complete the replacement of Dk by 
□ 9 in (|4T2"5]l , obtaining 

IMlLiCMo f ,) + ll v "llk s + ) + ll v "(^)lli 2 <~ [ V gU . S lU dV K 

(4.29) [ °' fo1 jA4 [o,io] 

+ ||Vu(0)||i 2 + ||x>7?n^||^ ( _ M ,+e sup \\Vu(t)\\ 2 L2 

' te[o,t ] 

To conclude the proof of (|4.23|) we estimate the integral term in the last inequality 
by i M . 1 1| □ 9 u||l£;* and use Cauchy-Schwarz. The last term on the right 

is eliminated by taking the suppremum in the resulting estimate over to G [0, t\] for 
arbitrary t\ > 0. 

b) We prove the estimate (|4.24[) for k = 1; the argument for larger k is identical. 
We begin by applying the estimate (|4.23j) to d{U. Commuting D g with &i we have 

O g dfU = d~ t U g u + 0(er _1_ )[(l - x P s)Q2U + XpsQiu] + 0(c(t))x ps Q2U 

where Q\ and Qi stand for second order operators with bounded coefficients. Here 
we have used (|4.8p for first derivatives of g away from the photon sphere, (|4.22[) for 
first derivatives of g near the photon sphere, and (|4.10|) for second order derivatives 
of g. We estimate the second term in \2 g dfU in LE^ and the third in L\L? X . This 
gives 

\\dM\hEl- ^AMlE^ 1 + \\ u \\lE^ + \\ a gU\\ LE ^i. 

Away from the event horizon the vector field di is timclikc, therefore, arguing as in 
the proof of Thcorcm l4.3f b). we can use an elliptic estimate to conclude that 

\\XouM\ le i k i < e\H\ LE ^ + IMIlb^ + \\ n gu\\ LE -^- 

On the other hand, near the event horizon we use Lemma 14.41 to obtain 

\\Xehu\\ H 2 < \\Xoutu\\ H 2 + \\XehVgU\\ H i + \\Wu(0)\\ H l 

Combining the last three estimates we obtain (|4.24[) for k = 1. □ 
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4.4. Conclusion. We can now prove Price's law for certain perturbations of the 
Kerr spacetimes: 

Theorem 4.6. Let g be a Lorentzian metric close to in the sense that it satisfies 
(|4.8[) . (|4.10l) and (|4.22[) . Let u solve (jl.ip with smooth, compactly supported initial 
data and V = 0. Then (fi~12j) holds. 

Proof. This is an obvious consequence of Theorems 11.51 14.31 and 14.51 □ 
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